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A homogeneous and isotropic cosmological model with a positive cosmological constant is con¬ 
sidered. The matter sector is given by a massless scalar field, which can be used as an internal 
time to deparametrize the theory. The idea is to study and compare the evolutions of a quantum 
and a classical probability distribution by performing a decomposition of both distributions in their 
corresponding moments. For the numerical analysis an initial peaked Gaussian state in the volume 
will be chosen. Furthermore, in order to check the robustness of certain results, as initial state both 
a slightly deformed Gaussian, as well as another completely different state, will also be studied. 
Differences and similarities between classical and quantum moments are pointed out. In particular, 
for a subset of moments classical and quantum evolutions are quite similar, but certain variables 
show remarkable differences. 
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I. INTRODUCTION 

During the last years an intensive effort is being made to construct effective theories (understood as a systematic 
framework that provides the classical equations of motion plus certain quantum corrections) for quantum cosmology 
[l|. The main motivation is to obtain testable results in scenarios where the complete knowledge of the underlying 
fundamental quantum gravity theory is not necessary, in such a way that quantum cosmology becomes an empirical 
science. In this respect, there are already several proposals and approaches in order to find quantum-gravity corrections 
to the anisotropy spectrum of the cosmic microwave background (see for instance M). 

In the particular case of loop quantum cosmology there are three key ingredients that one should consider in 
order to construct such an effective theory: holonomy corrections, inverse-triad corrections, and quantum-dynamical 
corrections. The origin of the first two corrections lies in the variables that are used in this specific theory: holonomies 
of connections and fluxes of spatial triads. Holonomy corrections appear in a process of regularization, and they take 
the form of higher powers of the connection that amend the classical Hamiltonian M- On the other hand, the 
inverse-triad corrections are produced because, in order to avoid infinities, the inverse of a given triad is replace d by 
the Poisson bracket between the corresponding triad and certain holonomy, which constitutes a classical identity [ll| . 
This procedure prevents, at the quantum level, the operator associated with the inverse of the triad from diverging, 
even when the triad itself tends to zero. The latter (quantum-dynamical) corrections are not specific to a loop 
quantization and arises due to the distributional character of quantum mechanics and the noncommutativity of the 
basic operators. In the present paper we will focus on these corrections, so let us analyze their origin in more detail. 

In the procedure of quantization each classical degree of freedom is replaced by an infinite set of quantum degrees of 
freedom, usually described by a probability distribution (wave function). Another way to parametrize these degrees 
of freedom is by decomposing the wave function in its infinite set of moments. These moments then appear in the 
classical Hamilton equations as quantum corrections. Nonetheless, this distributional character is not specific of the 
quantum theory. In fact, only ideally, classical mechanics is non-distributional. When the initial conditions are not 
known with infinite precision, this uncertainty can be described by an initial probability distribution, and it is then 
necessary to consider the evolution of such a distribution on the classical phase space. As explicitly shown in Bin, 
the evolution of a classical distribution can also be described in terms of its moments. And, interestingly, it turns out 
that the equations of motion for these classical moments can be obtained from the equations of motion for the quantum 
moments just by imposing a vanishing value of the Planck constant. This is a neat classical limit of a quantum theory, 
which implements the idea presented in 0 that the limit of a quantum theory is not a unique orbit on phase space, 
but an ensemble of classical orbits. Note that the mentioned classical and quantum probability distributions are 
defined on different spaces. Thus, their decomposition in moments allows us to compare their evolution. In addition, 
moments represent observable quantities that could be, in principle, experimentally measured. 
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Therefore, one can distinguish two different origins of quantum-dynamical corrections. On the one hand, the fact 
that an extensive (as opposed to a Dirac delta) distribution needs to be considered makes the presence of moments in 
the equations of motion unavoidable. Nonetheless, these kind of terms are also present in the evolution of a classical 
distribution and thus they are not genuinely quantum. For instance, as in the quantum case, they generically prevent 
the centroid of a classical distribution (the expectation value of the position and momentum) from following a classical 
trajectory on the phase space. On the other hand, the purely quantum or noncommutativity terms arise due to the 
noncommutativity of the basic operators. In the equations of motion they appear as a power series in Following 
the terminology of [l3j |. the first ones will be named distributional effects, whereas the latter ones purely quantum 
effects. It is important to stress that, if quantum cosmology is to become a testable theory, one needs to discriminate 
between purely quantum effects and effects that might appear just due to different technical or measurement errors, 
which would imply a classical probability distribution. 

Let us briefly mention that this formalism based on a decomposition of the wave function in terms of its moments 
was first presented in [l^ for the case of a particle on a potential. A similar formalism was derived in [l^, [l^ for 
generic Hamiltonians and on a canonical framework, but making use of a different ordering of the basic variables. 
Furthermore, this formalism has been adapted to the case when the dynamics is ruled by a Hamiltonian constraint, 
as opposed to a Hamiltonian function [i3- This is of particular importance in the context of general relativity, where 
the Hamiltonian is a linear combination of constraints. It has also been extensively applied to different models of 
quantum cosmology: isotropic models with a cosmological constant have been studied in whereas, in the 

context of a loop quantization, bounce scenarios have been analyzed in [^ . The problem of time has also been 
analyzed in [2l|, within this framework. Finally the classical counterpart of the formalism developed in EH, 
for the analysis of the evolution of a classical distribution in terms of its moments, was presented in [l3| and applied 
to the case of a particle on a potential in [^ . 

In this paper we will revisit, from the perspective of a classical probability distribution, the model studied in Ref. 
M for the evolution of quantum moments on a homogeneous and isotropic universe with a positive cosmological 
constant. The numerical analysis of the present model in terms of a wave function is presented in (^ . both for a 
geometrodynamical and a loop quantization. Here moments will be defined in terms of geometrodynamical variables. 

Our goal is to find differences and similarities between the quantum and classical (distributional) evolution of this 
cosmological model, given the same initial data for both cases, in order to study whether any moment has some 
distinctive or characteristic behavior under either the quantum or classical evolutions. As already commented above, 
the great advantage of the present formalism is that it considers the evolution of moments, which are observables, 
and no mention to an abstract mathematical object, as the wave function, has to be made. Therefore, choosing the 
same initial data for both classical and quantum sectors is straightforward. This kind of question is very difficult to 
be posed in terms of a wave function, since one needs to somehow define its classical probability distribution analog 
and choose the same initial data. (This can be done via, for instance, the Wigner transform but in general it is not 
positive definite and thus can not be strictly understood as a probability distribution.) Finally note that the present 
analysis is different from tying to define a classicalization of a given quantum system. This latter is usually addressed 
b^^pp ealing to the WKB limit or different decoherence processes that annihilate the quantum interference (see, e.g.. 

The rest of the paper is organized as follows. In Sec. |TT] the formalism for the evolution of classical and quantum 
probability distributions in terms of their corresponding moments is briefly summarized. Section fllll presents the spe¬ 
cific cosmological model under consideration. In Sec. llVI the results obtained with different numerical implementations 
are described. Finally, Sec. |V] discusses the main conclusions. 


II. GENERAL FORMALISM 

Let us assume a quantum mechanical system parametrized by the conjugate variables {q,p). The quantum moments 
associated to this system are given by 

G“'':=((p-p)“(g-g)Veyi, (1) 

where p := (p), q := (q) and the subscript Weyl stands for totally symmetric ordering. The order of a moment 
is defined as the sum between its two indices {a + b). Note that through this decomposition the wave function that 
describes the quantum state of the system gets replaced by its infinite set of moments G“’^, which only depend on 
time. Moments that correspond to a valid wave function must fulfill certain relations due to the Schwarz inequalities, 
see [T^ for a systematic derivation of such relations. The simplest, and probably most important, example is the 
Heisenberg uncertainty principle that, with this notation, takes the following form: 

(G1.1)2 <g2'0g°-2-^. 


(2) 
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The dynamical information of these moments is encoded on an effective Hamiltonian iJg, which is obtained by 
performing a Taylor expansion of the expectation value of the Hamiltonian operator around the centroid: 


OO OO 


1 na+& TT 

:= {H{q,p))weyi= {H{q q + q,p p + p))weyi = 


0=0 &=o 


a\h\ dp°-dq^ 


= H{q,p)+ E 


ci-\-b^‘2 


a\b\ dp°-dq^ 


( 3 ) 


In order to obtain the evolution equations for the moments as well as for the expectation values {q, p), it is 

enough to compute the Poisson brackets between each of these variables and the above Hamiltonian Hq. In this way 
one obtains the following infinite set of ordinary differential equations: 


dq 

dt 

dp 

dt 

dt 


dH{q,p) ^ 1 d-+^+^Hiq,p) , 

dp ^ a\b\ ’ 

^ o +&>2 ^ ^ 


dH{q,p) 

dq 


- E 

a+6>2 


1 d-+^+^Hiq,pl^^^, 


alb! dp°-dq^~^^ 


1 piC-\-d TT 


c+d >2 


c\d\ dp'^dq^ 


( 4 ) 

( 5 ) 

( 6 ) 


which is completely equivalent to the flow of states generated by the Schrodinger equation. Note that the first two 
equations are the usual Hamilton equations plus certain correction terms that depend on the moments. If all moments 
were vanishing, these terms would completely disappear. One of the consequences of these terms is that the centroid 
of a quantum distribution {q,p) does not follow a classical point orbit (the orbit obtained with an initial Dirac delta 
distribution, for which all moments vanish). Even so, this effect is not genuinely quantum since, as will be shown 
below, also happens for classical probability distributions. A closed formula for the Poisson brackets between any two 
moments, which has been left indicated in Eq. m, can be found in Ref. [Sill. For practical purposes, generically 
one needs to truncate the infinite system by introducing a cutoff by hand in order to be able to solve it. 

On the other hand, let us assume an ensemble on the classical phase space coordinatized by the conjugate variables 
{q, p). Such an ensemble will be described by a probability distribution function p(q, p, t ), which will obey the Liouville 
equation. This distribution function defines a natural expectation value operation on the classical phase space for any 
function f{q,p) in the following way. 


{f{q,p))c ■■= J dqdpf{q,p)p{q,p,t), (7) 

where the integration should be taken along the domain of the probability distribution. 

Following the same procedure as in the quantum case, the classical moments are then defined as 

G‘^’^:={ip-pnq-q)%, ( 8 ) 

where {q.,p) are the coordinates of the centroid of the distribution: q := {q)c and p := {p)c- [For notational simplicity, 
and for the moment being, the same notation {q.,p) will be used for the centroid of both classical and quantum 
distributions. Nonetheless, in Subsec. IHIBI specific notations will be introduced for each of them, to be used when 
the meaning is not clear from the context.] On the contrary to the quantum case, here all variables commute and, 
therefore, the ordering inside the expectation value is irrelevant. The Hamiltonian that describes the evolution of 
these classical variables is obtained, as in the quantum case, by expanding the expectation value of the Hamiltonian 
around the position of the centroid: 


Hc{q,P,C‘^’^) 


{H{q,p))c= H{q,p) + E 

Q-|-b^ 2 


1 

a\b\ dp°-dq^ 


(9) 


In order to get the evolution equations for the classical expectation values {q,p) and moments G“’^ one just needs 
to compute the Poisson brackets of different variables with this Hamiltonian. The flow generated by the Liouville 
equation for the probability distribution p{q,p,t) is then equivalent to the infinite set of equations obtained by this 
procedure. 
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It turns out that the difference between the quantum and classical equations are just terms that appear in the former 
equations multiplied by even powers of Planck constant and are missing in the classical equations. Such h factors 
are present in the quantum system due to the noncommutativity of the basic operators q and p. In fact, the classical 
equations of motion for the moments and expectation values can be obtained from their quantum counterparts by 
imposing a vanishing value of the Planck constant. Thus the classical limit of a quantum system, understood as 
h ^ 0, turns out to be very neat in this formalism. In particular, as can be seen, such a limit does not give a unique 
trajectory on the classical phase space but an ensemble of them. 

Due to the properties of the equations mentioned above, it is possible to classify the quantum effects depending 
on its origin. On the one hand, distributional effects arise because, due to the Heisenberg uncertainty relation, all 
quantum moments can not be vanishing. These effects are also present in the classical setting when considering the 
evolution of a spread probability distribution, for instance in the usual situation where the initial data are not known 
with infinite precision. On the other hand, the noncommutativity or purely quantum effects appear in the quantum 
equations as explicit h terms. As commented above, the origin of such terms lies in the noncommutativity of the 
basic operators. On the contrary to distributional effects, these are not present in the classical setting and are, thus, 
genuinely quantum. 

There are two classes of Hamiltonians that have very special properties regarding the classical and quantum evolution 
they generate [l3j |. On the one hand, the quantum equations derived from any harmonic Hamiltonian, which are at 
most quadratic on the basic variables, do not contain any h term. Therefore, the evolution that is generated by such 
a Hamiltonian both in the classical and quantum settings is exactly the same. On the other hand, the Hamiltonians 
that are linear in one of the basic variables, e.g. in g, generate the same evolution for the infinite set of variables 
(g,p, G"’°, G"’^) as for their classical counterparts (g,p, G^’^, G"’^) for all integer n. In particular, the cosmological 
model that will be considered on this paper is described by a Hamiltonian linear in q. In addition, when the 
cosmological constant is vanishing this Hamiltonian will turn out to be linear in both q and p. Thus, for that case, the 
Hamiltonian will be harmonic and the quantum and classical (distributional) evolutions it generates will be completely 
indistinguishable. 

Finally, let us comment that the stationary states can also be considered within the present formalism. The 
stationary states correspond to fixed points of the dynamical system under consideration and thus its corresponding 
moments can be obtained by solving the algebraic system of equations obtained by dropping all time derivatives from 
the equations of motion for (g,p, G“’^). In particular, following this procedure, the moments of the classical and 
quantum stationary states of the harmonic and the quartic oscillators were studied in [^ . 


III. APPLICATION TO A COSMOLOGICAL MODEL 


A. The classical cosmological model with positive cosmological constant 


Let us assume a cosmological model of a homogeneous, isotropic and spatially flat universe with a massless scalar 
field (j) as matter content and positive cosmological constant A. As usual in general relativity, this system is described 
by a Hamiltonian constraint, as opposed to a physical Hamiltonian. Nonetheless, it is possible to deparametrize 
the system and use the conjugate momentum of the scalar field which is a constant of motion, as a physical 
Hamiltonian. The Friedmann equation corresponding to this system reads 


a'V _ 47rGp| 


■A, 


( 10 ) 


where a is the scale factor and the prime stands for derivative with respect to the cosmic time. By choosing Newton 
constant as = 1 for convenience, it is straightforward to solve this equation for p^ and define our physical 
Hamiltonian as, 


H ■—prj, = a^\/|a'2 - Aa^l, (11) 

where the absolute value has been taken to extend this Hamiltonian to the region a'^ < Ka?. In this way, the scalar 
field (j) will play the role of time. Nevertheless, this Hamiltonian must still be written in terms of the canonical 
variables. Such variables are directly related to the scale factor of the universe as g = (1 — and p = —a^^a'. 

Then, the physical Hamiltonian takes the following form 


H = {\ — x)q\J\p‘^ — A[(l — 


(12) 
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The parameter x characterizes different possible cases of lattice refinement of an underlying discrete state, char¬ 
acteristic of a loop quantization [2^, [2^. A value around x = —1/2 seems to be favored by several independent 
phenomenological and stability analysis in the context of loop quantum cosmology [Si|. Nevertheless, in the 
(geometrodynamical) quantization studied here, with the usual canonical commutation relation, in principle this pa¬ 
rameter does not play any role. Thus, as was done in Ref. , X = —1/2 will be chosen since it turns out to be very 
convenient because it leaves the Hamiltonian as a linear function of the position: 


H = A|. 

The Hamilton equations take then the following form, 

q = 

P = 


3 - A) 




(13) 


(14) 

(15) 


where sg is the sign function and the dot stands for derivative with respect to 4>. The solution to these equations can 
be found analytically [l^ but here, in order to complement the discussion of that reference and to show in a more 
intuitive way the dynamics of this system in terms of the chosen variables, the corresponding phase space diagram is 
shown in Fig. [TJ With the chosen value of the parameter x, q is proportional to the volume of the universe a^, and one 
would then naturally define it as positive definite. Nonetheless, for illustrational purposes, both negative and positive 
values of the position q have been plotted. In fact, as can be seen in the diagram, the phase space is symmetric under 
a change of sign in q (as the whole system is symmetric under a change of sign of a). For positive q (as well as for 
negative q), the phase space contains three disjoint regions separated by the lines p = ±-\/A. The upper and lower 
sectors correspond to > A, whereas in the middle sector p^ < A. This latter region is the one that in principle is 
not allowed by the Friedmann equation (1101) (the momentum of the scalar field would need to be imaginary) but has 
been constructed extending the Hamiltonian by taking the absolute value inside the square root dm. 

The explicit equation for the orbits can be obtained analytically, 

(p2 - ^■)q^ = k, (16) 


with k an integration constant. Note that positive values of the constant k correspond to the upper and lower sectors, 
whereas negative values stand for orbits in the middle region. Between those solutions, there is also the degenerate 
solution <7 = 0, that corresponds to a universe with zero volume. 

In the upper region all orbits begin at the positive infinite of p and vanishing q, with a value (/> —>■ —oo of the 
affine parameter, and reach the asymptote p = \/A at an infinite value of q but at finite value of the affine parameter 
4> = </>div We will be interested in this region that corresponds to an expansion of the universe. The lower region is 
the time reversed of this latter one. There the system begins with an infinite volume and collapses. In fact, every 
orbit of this region is the analytic extension of the orbit with the same integration constant k of the upper region for 
values (j) G ((/div, c») of the affine parameter [^ . On the other hand, the middle region is nonphysical and corresponds 
to bouncing solutions that begin with infinite volume, collapse until a minimum value of g = ^J\k\/K with p = 0 and 
then expand again. Since we will be working just in the upper sector, from here on, the absolute values that appear 
in the Hamiltonian, as well as in the equations of motion, will be removed without loss of generality. 


B. Evolution equations for classical and quantum moments 

Following the procedure described in Sec. m one can construct the effective Hamiltonian Hq for this system as [h 


q ^_ q 

Hq = - A-f -VA^ 


n! 


gr„(p/VA)G"-° + n VAT„_i(p/VA)G”-1’^ 


(17) 


where the function 




(18) 
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FIG. 1: In this plot the phase space corresponding to the cosmological model under consideration is shown for the value A = 10. 
The dotted lines stand for the curves p = ±%/A, which can not be crossed by any orbit. These lines divide the phase space 
in three different regions. The upper region corresponds to expanding solutions. The lower region is the time reversed of the 
latter one and thus describes collapsing orbits. Finally, the intermediate region corresponds to bouncing solutions. 


has been defined. Making use of this Hamiltonian it is straightforward to obtain the infinite system of equations that 
rules the evolution of the expectation values {q,p) and quantum moments 


q = 


P = 


Qa,b ^ 


3 ^ A-"/2 


3 p 


(ZT„+i(p/VA)G"’°+ nVAT„(p/VA)G 


'n-1,1 


q __ q A —n/2 

- A- 


n\ 


^ nr ^ A-"/2 


n! 


qTnip/VX) {G“’^G”’°} +nyAr„_i(p/VA) {G“’^G”-^’l} 


(19) 

( 20 ) 
( 21 ) 


where the Poisson brackets between moments have been left indicated. (For explicit expressions of these brackets the 
reader is referred to [mil). The evolution equations for classical moments and expectation values {q,p,C°'’’’) are 
obtained from these previous ones just by imposing h = 0. Such h terms appear when computing the Poisson brackets 
between moments. 

Sometimes the meaning is not clear from the context thus, following the notation of the solution of the quantum 
system (HMD will be denoted as qq{(j)). On the other hand, the solution of the classical distributional system [that 
is, the one obtained from dlUSD by replacing all G“’^ by it corresponding G“’^ and imposing h — 0] will be denoted 
by gc(</')• Finally, the classical point trajectory [the solution to Eqs. dropping all moments] will be referred 

as qciass{4')- The same notation is used for the variable p. 

The system we are dealing with has a very particular Hamiltonian since it is linear on the position variable q. 
That is why only moments of the form G"’'^ and G"’^ appear in the effective Hamiltonian (HD, and in the second 
entry of the Poisson brackets of Eq. (ED, as well as in the equations for the expectation values dTM^Ol) . It can be 
shown that, for such a system, the infinite set of variables (g^, Pg, G"’°, G”’^) for all inte ger n obey a close system 
of equations, decoupled from the rest of the moments, which does not contain any h term [l3|. Therefore, given the 
same initial data, the expectation values g^, Pq, as well as all moments of the form G"’° and G"’^ follow exactly the 
same evolution as their classical distributional counterparts (gc, Pc, G”’°, G”’^). This is completely generic for any 
initial data. Therefore, for this kind of Hamiltonians gc('(') = g(j('(') and pd^f) = Pqi4’) for all times. However, this 
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does not mean that they follow the classical point trajectory, qciass{4’) 7^ since moment terms appear in Eqs. 

(fT5H^ni) . In the following subsections different initial data will be considered in order to compare the classical and 
quantum evolution of this cosmological model. And, thus, in order to find any difference, moments not contained in 
that subset will have to be checked. 

Finally, note that the A = 0 is a harmonic case since the Hamiltonian is linear in both position and momentum 
variables. This harmonic case has very special properties. In particular, all orders are decoupled, and there is no h 
in any of the evolution equation of the moments. Thus classical and quantum moments obey exactly the same set of 
equations. In fact, it is easy to find the analytic solution for all variables. 


9(0) 

P(0) 


9o exp 


^(0- 0o) 


Poexp - 

G a.b 
0 


^(0-0o) , 

^(&-a)(0-0o) 


( 22 ) 

(23) 

(24) 


{qo,Po, Gq’^) being the value of each function at tp = tpQ. 


IV. NUMERICAL IMPLEMENTATION 
A. Initial data 

In order to extract physical information from the infinite set of equations for quantum moments (UMII), as well 
as for its corresponding classical counterpart system, it is necessary to resort to numerical methods. In addition, for 
practical purposes, one needs to consider a cutoff, that is, a maximum number N for which all moments of an order 
greater than N are dropped. This fact shows that this method of moments is well suited for peaked states, when 
high-order moments are negligible. Nonetheless, the usual tendency of the dynamics is to spread out the state so 
that, from certain point on, this method will not give trustable results. There are several control methods to know 
when this happens. On the one hand, one should study the convergence of the solution with the cutoff. This is 
done by solving the system of equations with different cutoffs and checking that the difference between solutions with 
consecutive cutoffs tends to zero. On the other hand, the conservation of constants of motion should also be taken 
into account. In the present model the Hamiltonian Hq (and He for the classical system) is conserved. Finally, the 
high-order inequalities obtained in should also be fulfilled during the whole evolution. All these control methods 
have been used to test the numerical implementations that will be presented below. 

The goal of the present paper is to compare the evolution of the quantum and classical (distributional) systems 
within the context of the formalism of moments described above. With that purpose, similar numerical evolutions to 
those that were performed in Ref. will be performed but in this case not only for quantum moments, also for their 
classical counterparts. In that reference a Gaussian state was chosen as initial state. The moments corresponding to 
such a state read, 


Qa,b 


2-{a+b)f^a^b-a ^ ^ 

I 2 j'U J- 


0 otherwise, 


(25) 


a being the width of the Gaussian. As can be seen, this state is of a very special kind, as many of its moments vanish. 
In fact, in order to check whether the results obtained depend on the fact that these moments are vanishing; apart 
from the Gaussian initial state, here another two initial states will also be considered. In other words, three evolutions 
will be performed. The first one with a Gaussian state for both classical and quantum moments. Whereas the initial 
state for the second evolution will be given by evolving the Gaussian state with the quantum equations during a 
short period of time (O.l^div). At this point all moments will have been excited and the resulting state will be a 
slightly deformed Gaussian, which will be used as initial state for both classical and quantum systems. Obviously, 
both evolutions will give the same result for the quantum moments but not for the classical ones. Finally, the initial 
data for the third evolution will be given by all classical and quantum moments taking the nonvanishing value 

G“’'' = alblh"^, 


(26) 









which is allowed by all high-order uncertainty relations obtained in [l3l| . This is a distribution quite different from 
the Gaussian (051) . which will be used to check the generality of the results obtained in the previous two cases. 

Regarding more technical issues, the numerical evolutions have been performed for all cutoffs from = 2 to 

= 10 in order to verify the convergence of the method with the cutoff order. Furthermore, different values of the 
cosmological constant have been considered: small (A = 1), intermedium (A = 10"^), and large (A = 9 x 10^). The 
initial conditions for the position and the momentum have been chosen as p{0) = 10®, q{0) = 1 to ensure that the 
state is on the upper region of the phase space shown in Fig. [I]and corresponds to an expanding solution. In addition, 
the Gaussian width has been taken as cr = so that the expression for a moment G“’® is completely symmetric on 

its indices a and b. In this way, initially = G®’“ and the fluctuation (as well as higher-order moments) of the 

position are equal to those of the momentum. Finally, the numerical value of the Planck constant has been chosen as 
A = 10-2. 

In summary, the main difference with respect to the Gaussian state used in Ref. [l^ is that its width was chosen as 
a = 10-2, that is, smaller than in our present case with a = 10-^. In addition, in that reference, h was chosen to be 
equal to one in order to make the backreaction effects more clearly visible. However, qualitatively the results do not 
change with such modifications. Here, in order to analyze the purely quantum effects, which appear in the quantum 
equations of motion as a power series in h?, it is necessary to consider a smaller Planck constant, otherwise all terms 
of the form would be of the same order. 

In the next subsection we will comment, mainly qualitatively, the behavior of different variables throughout evolution 
for the Gaussian initial state (1251) for both classical and quantum moments; whereas in Subsec. IIV Cl the differences 
between classical and quantum moments will be analyzed quantitatively in detail. 


B. Description of the evolution 

Regarding the expectation values, as has already been commented in the preceding section, due to the linearity of 
the Hamiltonian they have exactly the same evolution both for the classical (distributional) case as for the quantum 
case, that is, = qc and Pq = Pc at all times. Nonetheless, the backreaction on the classical orbits is not vanishing due 
to the presence of moments in Eqs. (muni)- In other words, the centroid of a (classical or quantum) distribution will 
not follow a classical point trajectory on phase space. This backreaction can be measured by defining the differences 
Sq := qq — qdass and Sp := Pq—pdass- [Note that in principle qq{t) should be the solution of Eq. (1191) by considering the 
infinite set of equations. We obtain an approximate version of this solution by imposing the cutoff N = 10.] Regarding 
this backreaction on the classical orbits, though less severe due to the smaller numerical value of h, we obtain similar 
results as in Ref. [l^ . In particular the deviation from the classical orbit is greater as we move to larger values of the 
cosmological constant. More precisely, measured at 0.8(/)div, the relative change in volume is Sq/q ~ IQ-^ for the large 
cosmological constant case. In the other two cases, due to its smallness, this effect is mixed with numerical error and 
it turns out to be difficult to measure but, in any case, we have that Sq/q ^ IQ-^® throughout evolution. In all cases 
there is enhancement of the divergence in the sense that the position qq (or q^) corresponding to the centroid of a 
probability distribution approaches the divergence faster than the classical point orbit qdass- On the other hand, the 
deviation (in absolute value) of the momentum p is not that pronounced and for all cases, during the whole evolution 
until the mentioned time, Sp/p ^ IQ-^®. Thus, in this sense, the case with larger cosmological constant A = 9 x 10^ 
is the one where the backreaction effects are more relevant. The evolution of the volume q is represented in Fig. [2] 
for different values of the cosmological constant. 

Regarding the evolution of different moments, as already commented above, due to the linearity of the Hamiltonian, 
pure fluctuations of the momentum G”’° and moments of the form G"’^ coincide with their classical counterparts: 
Qn,o _ (jn,o (^n,i _ fjn.i times. For the rest of the moments it can be asserted that, the behavior 

of classical moments G“’® is qualitatively similar to their quantum counterparts G“’® except for those of the form 
G°^2 "+i, jn the next subsection we will analyze in more detail the differences between both classical and quantum 
moments, and particularly the corresponding to these pure-odd fluctuations of the position G®’2"+^. For the rest of 
this subsection, and unless explicitly stated, all that is said about the classical moments G“’® applies equally to the 
quantum ones G“’®. In summary, the generic behavior of the moments that will be commented below can be seen in 
Fig. [3] for the particular case of two moments: G°’2 and G^’®. 

For a Gaussian state all moments with an odd index are vanishing. These moments get excited as soon as the 
evolution begins because their time derivative is nonzero. After that excitation generally moments |G“’®(()))| with a > b, 
including those that have been excited from an initial vanishing value, behave as an (approximately exponentially) 
decreasing function; whereas moments with a < b increase in absolute value. This behavior is inherited from the 
solution ((331) for the harmonic case A = 0. [Nevertheless, even if it is very convenient to have such a picture, several 
moments break this rule (usually moments that are supposed to be decreasing turning out to be increasing with time), 
specially as the value of the cosmological constant is larger.] In fact, this general tendency can be intuitively expected 
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log,„9 



FIG. 2: In this figure the evolution of the volume q is plot¬ 
ted for three different values of the cosmological constant on 
a logarithmic scale. The green (continuous) line stands for 
A = 1, the blue (dot-dashed) line corresponds to A = 10^, 
whereas the black (dotted) line represents the solution with 
A = 9 X 10^. Note that for the same value of 
q{(j>/(j>div) is larger, the smaller the value of the cosmological 
constant. 


log,„C»Mog,„|C''’| 



FIG. 3: In this plot the evolution of the moments and 
is shown in a logarithmic scale for three different values 
of the cosmological coustant and for an initial Gaussian state. 
Gontinuous lines correspond to C^’®, which is an initially van¬ 
ishing moment, whereas dashed lines stand for C°’^. The color 
(thickness) of the lines represent different values of A: green 
(thickest) for A = 1, blue (intermedium) for A = lO'*, and black 
(thinnest) for A = 9 x 10^. Note that the larger the value of 
A, the larger the slope of the moment in terms of 


by looking at the phase space of the system shown in Fig. [TJ Let us assume that initially we have a homogeneous 
and compact probability distribution with the shape of a circle centered at small qq and large pq. When the evolution 
begins, each point of that circle follows its corresponding orbit and the circle will get deformed. In particular, as these 
points approach the divergence at p = \/A, the initial circle will be enlarged in the horizontal (g) direction whereas it 
will be contracted in the vertical (p) direction. Thus, the state will be more spread in the q direction but more peaked 
in the p direction. This is exactly what we have described in terms of moments: moments (7“’^ with more weight in 
the p direction (a > h) will be decreasing, and the others increasing, in absolute value. 

On the other hand, regarding the sign of the moments, during the initial stages of the evolution we observe that, 


^even,even ^ q ^even,odd ^ q 

^odd,even ^ q ^odd,odd ^ q 


(27) 

(28) 


The only exception to these rules are the pure-odd fluctuations of the momentum which are positive. The 

same applies to quantum moments, except for those of the form that, as will be explained below, are initially 

excited to a negative value and then change sign during the evolution. Note that the first of the inequalities, which 
states that moments with both even indices must be positive, is implied by the very definition of the moments and 
must be obeyed at all times for any system [l3l| . 

Even if the commented qualitative features for the moments are independent of the value of the cosmological 
constant, there is indeed some dependence. In particular, the larger A the lower the slope, in terms of 4)/4>div, 
of a given moment is. [Note, however, that t/fdiv also depends on the cosmological constant since ^div := + 

2/3tanh~^(y^Pg — A/po)]. In addition, the excitation value of the initially vanishing moments increase, in absolute 
value, as we consider larger values of the cosmological constant. However, since the slope of the moments is larger 
for small A, at the end stages of the evolution an increasing (decreasing) moment has usually a larger (smaller) 

value for smaller values of A, as can be seen in the example of Fig. [3] Even so, the relative fluctuations 0°“’^/{q°‘p^) 
are generally larger the larger the value of the cosmological constant. 


C. Quantitative comparison between classical and quantnm moments 

Following the notation of [l^, we introduce the following operators, and ^ 2 , to quantitatively measure the 
difference between the classical and quantum evolution of the present system: 


6iq{(l)) := qci4>) - qciass{4>), 

^2q{(t>) ■= qqif/}) - qM)- 


(29) 

(30) 
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Note that the action of i5i on a given moment is not defined since there are no moments in the classical 

point orbit; whereas 620 °"’^ := — G“’^. The first operator can be understood as a measure of distributional 

effects, whereas the second one (5i measures the strength of purely quantum effects. 

Due to the properties of linear Hamiltonians, in this case ^29 and 52 P are vanishing. Therefore, the departure 
(Jig, 6 ip) of the centroid from the classical (point) trajectory on the phase space is uniquely due to distributional effects, 
which are indeed exactly reproduced by a similar distribution evolving on the classical phase space. Furthermore, 
J 2 G"’° and J 2 G”’^ are also vanishing for any integer n given any initial data. Therefore, in order to check the relevance 
of the purely quantum effects of this system, it is necessary to consider other kind of moments. 

In order to check the robustness of certain results that will be commented below, three different sets of initial data 
have been considered: a Gaussian state with moments given by (1^ . a slightly deformed Gaussian (by evolving the 
previous Gaussian with the quantum equations during O.l^div)) and a state with all nonvanishing moments of the 
form (1^ . As one would expect, for the deformed Gaussian case, the classical moments tend to show less divergence 
from their quantum counterparts as in the Gaussian case. 

The analysis of the results of such numerical implementations has shown several interesting features, which are listed 
below from i/ to iv/. Note that, in order to remove possible spurious effect of the cutoff (iV = 10), only moments up 
to order seven will be considered for the present discussion since higher-order moments are the most sensitive ones 
to the effect of the cutoff. Unless otherwise explicitly stated, the following results correspond to the Gaussian initial 
state. In particular, we will refer to the implementations of the other two initial data at the final part of item in/. 

i/ One of the important results is that all moments G“’^, except those of the form even if quantitatively 

different, have the same qualitative behavior as their classical counterparts G“’^. Thus purely quantum effects act, as 
one would expect, as small perturbations by slightly deforming the numerical values of different variables but not, in 
general, its qualitative behavior. 

a/ In absolute terms, |J 2 G“’^| is largest for moments of the form G°’” and G^’^". This absolute difference tends to 
increase with time and depends on the value of the cosmological constant. More precisely, for the case with A = 1 
at a time </) = O.750div, the highest difference corresponds to the moment G°’® with J 2 G°’® « 6 x 10"^. Such a large 
difference is only measured for this moment, which is also the largest of all considered ones [G°’®(0.75(/)div) ~ 4 x 10 ^^]; 
the next highest value being J 2 G°’^ ~ —7, whereas the rest of the moments have an upper bound |J 2 G“’^| < 0.1. On 
the other hand, for the intermediate value of the cosmological constant A = 10^ at the mentioned time </> = 0.75(/)div) 
all absolute differences are |J 2 G“’^| < 6 x 10“'^, whereas for large cosmological constant value A = 9 x 10^ this upper 
bound is much lower: |J 2 G“’^| < 2 x 10 ”^. 

iii/ In order to measure the strength of purely quantum effects in relative terms, we define the relative difference 
Ja,& := J 2 G“’^/G“’^ and the ratio Va.b '■= G“’^/G“’^, which is probably more intuitive to understand. These two objects 
are obviously related as ra.b = 1 — 

For the Gaussian initial state, in relative terms, we find that moments G“’^ with a + 6 an odd number, that is, 
with an even and an odd index, suffer the largest departure from their classical counterparts G“’^. [This statement, 
of course, exclude those of the form G"’° and G"^.] For the remaining moments, even if generically increasing (in 
absolute value) with time, so that the largest values are measured at the end of the evolution, their relative changes 
are bounded as |Ja,b| < 4 x 10“^° for the cases A = 1,10'*; whereas for the case of a large cosmological constant 
A = 9 X lO*’, jJa.bl < 2 X 10-6. 

Let us then focus on these moments with largest departure, that is, moments G“’^ with a + 6 an odd number. The 
evolution of their ratios Vafi is plotted in Fig. 0] for different values of the cosmological constant and for both the 
Gaussian (left column) and the deformed Gaussian (right column) initial states. For the Gaussian initial state their 
behavior can be separated in the following three different kinds: 

1 . Moments of the form As already commented above, a given moment of the form general does 

not follow, even qualitatively, the same evolution as its quantum counterpart ( 70 . 271 - 1-1 _ come back to 

these latter in the point iv/. 

2. Moments of the form (^odd.even j^emai-kably, as can be seen in Fig. |4l the relative change of these moments 
are constants of motion regardless the value of the cosmological constant. The exception to this rule is G^’^ 
which, for all values of the cosmological constant, is initially excited to a value much higher than its classical 


counterpart G®’^. Thus their ratio rb ,2 is very small but increases slowly with time, 
of the moments of this type take approximately the following values: 

The ratios r^^b for the rest 

ri,2 ~ ?’i,4 ~ ri^G Ki 2/3, 

(31) 

r3.2~ 1 / 2 , 

(32) 

r3.4~ 2/5. 

(33) 


All these moments are negative and increasing (in absolute value) throughout evolution. 
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FIG. 4: In these plots the evolution of the ratios ra,b = Q°-,b^Qa,b corresponding to quantum moments with highest deviation 
from their classical analogs, excluding those of the form jg shown. The left column corresponds to the initial Gaussian 

state, whereas the right column stands for the initial deformed Gaussian case. The upper plots correspond to A = 1, the 
medium ones to A = 10“*, and the lower ones to A = 9 x 10^. The (continuous) black line represents the three ratios ri, 2 , ri, 4 , 
and ri,6. All of them appear superposed since they follow very similar trajectories. The ratios r 3,2 and r 3,4 correspond to the 
red (dotted) and green (long-dashed) lines respectively. Finally, r 2,3 appears represented by the gray (continuous) line, r 2,5 by 
the black (dashed) line, and r4,3 by the orange (dot-dashed) line. Note that the late-time behavior is the same for both initial 
states. 


3. Moments of the form G'even.odd^ These moments, which, up to seventh order, are just three (G^’^, G^’®, and 
G^’^), are positive and increasing (in absolute value) during evolution. They are initially excited to a higher 
value than their classical counterparts, thus Va^b > 1- Remarkably, this initial value of the ratio is independent 
of the value of the cosmological constant. Nonetheless, and contrary to the previous case, the relative difference 
between classical and quantum moments does depend on time. More specifically, ratios ra,b tend to decrease 
with time, which means that the classical moments increase faster than the quantum ones. Furthermore, this 
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decrease is faster the larger the value of the cosmological constant. Interestingly, r 2,5 and r 2,3 seem to tend 
to the same asymptote, as is more clearly visible in the last plot of the left column of Fig. 0] For the case 
with largest cosmological constant (A = 9 x 10^), the final values of all these ratios are lower than one, thus 
the corresponding quantum moments G“’^ end up being smaller than their classical counterparts. On the other 
hand, for the other two cases (A = 1 and A = 10"^), the final values of these ratios are still larger than one. 

Note that all moments which have been analyzed in the last three points were initially vanishing. Hence, we can 
divide the analysis in two stages: the initial excitation and subsequent evolution. All these moments are excited to 
different values than their classical counterparts. Therefore, there is a clear purely quantum effect acting on this 
excitation mechanism (that is no other than the h terms present in the equations of motion). Nevertheless, once this 
mechanism acts we observe different behaviors. On the one hand, those moments that have been initially excited to a 
smaller value than their classical counterpart {ra,b < 1) keep their relative difference constant throughout evolution, 
that is, they evolve in a very similar way as their classical counterparts. This means that, for this kind of moments, 
purely quantum effects act continuously during the whole evolution such that is always proportional to G“’^, 
with a constant proportionality coefficient. On the other hand, the evolution of the quantum moments that have 
been initially excited to a larger value than their classical counterparts {ra,b > 1) is slowed down by purely quantum 
effects. In such a way that the corresponding classical moments increase faster (in absolute value). Therefore, the 
net effect of the purely quantum terms is contrary in the initial excitation of the quantum moment, when its value is 
risen with respect to its classical counterpart, and during the subsequent evolution, when it is slowed down. 

These results are quite unexpected, especially the constant behavior of the ratio between certain moments. In fact, 
this was the main motivation to consider other sets of initial data. Regarding the initial deformed Gaussian, the 
results obtained are shown in the right-hand column of Fig. |4l It is clearly seen that initially the ratios Va^b behave 
in a different way as in the previous case, in particular none of them is kept constant. Nonetheless, after a period 
of time, the classical moments follow the same tendency they had for the Gaussian case and the final part of the 
evolution regarding the different ratios have exactly the same shape as before. 

In addition, just to test that this behavior is not completely generic, as would be expected, for this particular 
issue another third set of initial data has been considered, given by Eq. (E5)) . In this case, the behavior described 
above for the ratios f, of moments of the form codd.even Qeven.odd disappear and the same tendency as the 
remaining moments is shown; that is, all ratios are close to one initially and they depart from one as the evolution 
advances. Nevertheless, this departure is small and all relative differences are bounded as liJa.bl < 5 x 10“® for A = 1, 
< 4 X 10“”^ for A = 10^, whereas for the the large cosmological constant A = 9 x 10^ case |(5o,&| < 5 x 10“^. 

Thus, we conclude that the behavior found for the Gaussian case is robust under small deformation of the initial 
data, but not completely generic. This situation, nonetheless, might be important if we know that the initial state is 
an approximately semiclassical Gaussian state. 

In any case, finding an analytical explanation of this surprising behavior is very difficult since the system under 
consideration is a nonlinear and highly coupled set of differential equations. As an example, in the appendix the 
purely quantum terms that appear in the equations of motion for the moments G^’^ and G^’^ are shown. Even if 
the terms are quite similar, in the sense that the same moments (except one) appear in both equations, the behavior 
of their corresponding ratios r 2,2 and r 3 ^ 2 ) which somehow quantifies the effect of the presence of these terms, are 
completely different. 

iv/ Finally, let us analyze the behavior of moments of the form G°’°‘^'^. As has been commented above, these are 
the only moments that follow qualitatively different trajectories in the classical and quantum settings. Let us be 
more specific. For the Gaussian initial state in the classical case all moments G°’°'^‘^, which are initially vanishing, 
are excited to certain positive value and increase through evolution. On the contrary, their quantum counterparts 
GO’Odd are initially excited to a negative value (much smaller in absolute value than their classical counterparts) and 
their evolution is decreasing (that is, increasing in module), which makes their corresponding ratio Va^b negative but 
close to zero. This decreasing behavior follows during the whole evolution for the cases of small and intermedium 
cosmological constant cases, which means that quantum moments increase faster (in absolute value) than their classical 
counterparts. Nonetheless, for the large cosmological constant case A = 9 x lO"^, at certain point these moments get 
a minimum and increase afterwards. Remarkably at the same time, around (p = 0.346(;iidiv, ah moments Gd’°dd cross 
zero and follow their growing behavior (see Fig. [S]). Finally, each quantum moment Gd’^"+^ tends asymptotically 
to its classical counterpart Gd’^"+^, as can be seen in Fig. [SI which represents the evolution of the ratio ro ,3 as an 
example of moments of this kind. For the deformed Gaussian moments are initially negative. Nevertheless all classical 
moments cross zero and tend to the above described evolution for the initial Gaussian case. 
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FIG. 5: The initial stages of the evolution of the moments 
(red dot-dashed line), (red continuous line) and 
their classical counterparts G°’^ (black dot-dashed line), 
G°’® (black continuous line) for the particular case with 
A = 9 X 10^. 


0.3 



FIG. 6: The ratio ro^s is shown for different values of the 
cosmological constant. The blue (thickest) line corresponds 
to A = 9 X 10^, the brown to A = 10"^, and the red 
(thinnest) line to A = 1. The dotted black line just rep¬ 
resents the asymptote ro, 3 =l. Note that, for small and in¬ 
termedium values of the cosmological constant, ro ,3 is neg¬ 
ative throughout evolution due to the fact that G°’^ is neg¬ 
ative, whereas G°’^ is positive. On the contrary, in the case 
A = 9 X 10^, G°’® flips sign, which produces the vertical 
line around (p = 0.346())div, and the ratio ro ,3 tends to one 
asymptotically. 


V. CONCLUSIONS 

In this paper the formalism presented in [T^ to study the evolution of classical and quantum probability distribu¬ 
tions, by performing a decomposition on its corresponding moments, has been applied to a particular cosmological 
model. The goal was to study similarities and differences between the dynamics of classical and quantum moments, 
in order to find physical consequences of purely quantum terms, which are present in the equations of motion for 
quantum moments. 

The cosmological model under consideration is a homogeneous and isotropic universe with a massless scalar field 
and a positive cosmological constant. The evolution of such a model is ruled by a Hamiltonian constraint, which 
can be deparametrized by using the scalar field as internal time and, thus, its conjugate momentum as the physical 
Hamiltonian. After introducing an appropriate pair of conjugate variables, that represent the volume and the Hubble 
factor of the universe, this physical Hamiltonian turns out to be linear in the volume q. 

In fact, it can be shown that for any linear Hamiltonian in q the evolution of the expectation values {q,p), as well 
as of the moments of the form (G"’°,G"’^) for any integer n, coincides both in the quantum and classical settings. 
In other words, given the same initial conditions, the centroid and the mentioned moments of a classical distribution 
evolve in exactly the same way as the corresponding variables of a quantum distribution. Therefore, in order to find 
physical consequences of the purely quantum effects, the dynamics of moments not contained in that subset needs to 
be analyzed. 

Nonetheless, for such a purpose, it is necessary to resort to numerical methods due to the complexity of the equations 
of motion. In particular, three different initial states have been considered for both classical and quantum settings. 
First a Gaussian in the volume has been chosen. Second a slightly deformed Gaussian state has been constructed 
by evolving the previous Gaussian state during a short period of time with the quantum evolution equations. The 
obtained state has then be used as initial condition for both classical and quantum moments. Finally, the third set of 
initial data is given by Eq. (j26l) and represents a state which is not close to a Gaussian. This state has been mainly 
used to check that the behavior of the ratios shown in Fig. |T]is not completely generic. 

The main result is that, in relative terms, moments of the form G“’^ with a -I- & an odd number, are the ones that 
show most divergence from their classical analogs. For the initial Gaussian state, all these moments are vanishing. 
Since their time derivative is nonzero they are immediately excited as soon as the evolution begins. This excitation 
value turns out to be quite different for a given classical moment and its quantum counterpart. After this excitation, 
two different behaviors have been observed. Quantum moments that have been excited to a smaller value than their 
classical counterparts (va^b < 1) evolve in a very similar way as their classical analog, keeping the ratio ra,b constant. 
On the other hand, evolution of the quantum moments that have been initially excited to a larger value than their 
classical counterparts {ra^b > 1) is slowed down by purely quantum effects. In this way, their corresponding classical 
moments grow faster in absolute value. In summary, the effect of the purely quantum terms is contrary initially and 










14 


during evolution. While initially the value of the quantum moment is risen with respect to its classical analog, during 
the subsequent evolution its growth is slowed down. 

In addition, for the initial deformed Gaussian state, even if initially different, the ratios ra,b tend to the same 
behavior as explained above for the Gaussian case. On the contrary, the third set of initial data (1251) does not follow 
this tendency. In this latter case, all ratios depart only slightly from one. Therefore, we conclude that this behavior 
of the different ratios is robust under small deformations of the initial Gaussian state but, certainly, not completely 
generic. 
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Appendix A: An example of purely quantum terms 

In this appendix, the purely quantum terms for the equations of the moments and G^’^ are shown. These are 
the terms with an explicit h factor that appear in the right-hand side of their corresponding equations of motion: 

G2.2 ^ 17/2 |48p (/ - G2’° - 24 (a + 4/) (/ _ g^,o ^ ^ 4 ^ 2 ^ (^2 _ ^^4 ^ 4.0 

64 (jp — A) ' 

- 30 (A^ -k 8 / -k I2A/) (p2 _ A)^ GAO -k 42p (SA^ -k -k 20Ap^) {p^ - A)^ G®’° 

- 7 (5A3 -k 64p® -k 240Ap4 -k I20A2p2) (^2 _ ^ 7.0 (- 35^3 54^6 330Ap4 2S0A^p‘^) G®’°| -k ... 

G3.2 ^ / - 16 _ a)"g2-o + 48p - Af G^’O - 24 (A + 4p") ip^ - A? G^o 

128 (p 2 - I ' 

-k 40p (3A -k 4p2) (p 2 _ A)'^ G'^’° - 30 (A^ -k 8 p^ -k l2Ap^) {p^ - A)^ G®’° -k 42p (SA^ -k 8 p^ -k 20Ap^) {p^ - A)^ G^’° 

- 7 (5A3 -k 64p® -k 240Ap4 -k 120A^p'^) {p^ - A) G®’° -k 9p (SSA^ -k 64^^ -k 336Ap‘^ + 280A^p"^) G®’°| -k ..., 

where dots stand for terms with no h factors. As can be appreciated, the purely quantum terms that appear in these 
two equations are pretty similar, in the sense that the same moments appear (except one) with similar polynomial 
coefficients depending on the expectation values q and p, as well as on the cosmological constant. Nonetheless, the 
effect of such terms are very different on each case. On the one hand, G^’^ is quite similar as G^’^ during the whole 
evolution, keeping r 2,2 around one but increasing with time. On the other hand, for an initial Gaussian state, 73^2 
takes a constant value around 1/2, independently of the value of the cosmological constant. The behavior of the latter 
ratio r 3_2 can also be seen in Fig. Ufor the evolved Gaussian initial state. In that case, it begins with a value of one, 
since both quantum and classical moments are initially equal and decreases until reaching again the constant value 
of 1/2. 
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